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THE LINEARIZED CHARACTERISTICS METHOD AND ITS APPLICATION TO PRACTICAL 

NONLINEAR SUPERSONIC PROBLEMS 1 

By Antonio Terri 


SUMMARY 

The method of characteristics has been linearized by assuming 
that the flow field can be represented as a basic flow field determined 
by nonlinearized methods and a linearized superposed flow field 
that accounts for small changes of boundary conditions. The 
method has been applied to two-dimensional rotational flow 
where the basic flow is potential flow and to axially symmetric 
problems where conical flows have been used as the basic flows . 
In both cases the method allows the determination of the flow 
field to be simplified and the numerical work to be reduced to a 
few calculations. The calculation of axially symmetric flow 
can be simplified if tabulated values of some coefficients of the 
conical flow are obtained. The method has also been applied to 
slender bodies without symmetry and to some three-dimensional 
wing problems where two-dimensional flow can be used as 
the basic flow. Both problems were unsolved before in the 
approximation of nonlinearized flow. 

INTRODUCTION 

The use of the method of characteristics for the solution, of 
supersonic-flow problems requires numerical procedures 
which are lengthy and involved and which must be repeated 
for each set of boundary conditions. The method has received 
general practical application only for two-dimensional or 
axially symmetrical problems in steady flow and one- 
dimensional or quasi-one-dimensional nonsteady flow, and 
only very few eases of general three-dimensional flow have 
yet been investigated. 

Man}” problems have been investigated at present by 
means of the linearized theory in which the disturbance- 
velocity components (defined as the difference between the 
local and the free-stream components of the velocity) are 
considered and are small, so that terms of second order or 
higher can be neglected. In the present report a simplifica- 
tion is introduced in the equations of motion based on the 
assumption that one of the velocity components or the 
variation of the velocity components as a function of a given 
parameter can be considered small, so that terms of second 
or higher order in the quantities considered small can be 
neglected. When one of the velocity components is assumed 
to be small, the other two velocity components can be ex- 
pressed in two parts, one of which is large and is a function 
only of two coordinate positions, and the other of which 


is small, of the same order as* the third velocity component, 
and is a function of all three coordinates. 

If the variations of velocity components as functions of a 
given parameter are considered small, all three velocity 
components can be expressed in two parts. One, large, is 
independent of the parameter considered, and the other, 
small, is a funetipn of the parameter considered. When the 
velocity components are substituted into the differential 
equations, the equations can be divided into two parts, and 
the differential equations containing the velocity components 
considered small become linear; therefore, superposition of 
solutions .is possible. With this assumption the Sow field 
can be represented for any condition by the superposition 
on a nonlinear basic flowfield of a linearized flow perturbation. 
The flow field, which represents the variation of the basic 
flow due to the changes of the geometrical parameter con- 
sidered, changes linearly with the parameter. Because of 
the simplification, the superposed flow field is defined by 
differential equations of hyperbolic type which have char- 
acteristic surfaces equal to the characteristic surfaces of the 
basic flow field and known coefficient; therefore, the super- 
posed flow field can be obtained directly without the iteration 
process along the characteristic net of the basic flow. 

A particular application of the linearized characteristics 
method has been discussed in references 1, 2, and 3 in which 
bodies of revolution at small angles of attack have been 
considered. In the present report the basic concept of the 
linearization is discussed and examples of application to two- 
dimensional rotational flow, to conical flow, to axially sym- 
metric flow, and to some general three-dimensional problems 
are discussed. From these examples, other applications of 
the same method to supersonic steady- or nonsteady-flow 
problems can be visualized. For example, the method can 
also be applied to the determination of the flow field in super- 
sonic compressors or turbines having supersonic relative 
velocity inside the passage. In this case, the two-dimensional 
flow of the cascade of the blades at each radial station 
or the axially symmetric flow can be assumed as the basic 
flow. In the first case the radial component of the velocity 
must be assumed to be small, whereas, in the second case, 
the tangential component of the velocity must be assumed 
to be small. 


■ • 

1 Supersedes MAC A TN 25L5 t “The Linearized Characteristics Method and Its Application to Practical Nonlinear Supersonic Problems” by Antonio Ferrt, 1651. 
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SYMBOLS 


a, b, c, d, e 


A n 

B and C 

C, 

M 

R 

S 

U; V, w 


tv, V n , w 
v T , Vn 

V 

P 

y 

Xi.a 

d 


<P 

v> i" 


coefficients of the variables in the characteristic 
equations (defined case by case) 
speed of sound 

coefficient ( j~^ (ttoWn+Po p,)^ 

coefficients defined by equations (20) and (25) 
pressure coefficient 
Mach number 

gas constant or radius of hodograph diagram 
entropy 

velocity components along the x-, y-, and 2 -axes 
in Cartesian coordinates or along x- and y - axes 
and perpendicular to the meridian iry-plane 
in cylindrical coordinates 
velocity components in polar coordinates 
tangential and normal velocity components in 
front of the shock 
intensity of the velocity vector 
Mach angle 
ratio of specific heats 

inclination of the characteristic lines in the plane 
2 =Constant or 0=Constant 
coordinate of the meridian zy-plane in cylindrical 
coordinates or of the r^-plane in polar 
coordinates 
polar coordinates 

inclination of the velocity vector with respect to 
the a;- ax is 

components of the rotation along the x-, y-, and 
2 -axes 


Subscripts: 

0 properties of the basic flow 

1 . . . n properties of the superposed linearized flow fields 


THE EQUATIONS OF THE LINEARIZED CHARACTERISTIC 

SYSTEM 


Consider, for example, a flow field defined by a velocity 
vector V(u,v,w) the components of which can be expressed in 
the form 


N “) 


A r 


t 


y 


y 

i 


( 1 ) 


Assume that «o and are functions only of x and y; then, by 
neglecting terms of the order of a s 2 , equation (2) becomes 


bS a 2 {bv 0 bu 0 \ „ {be, bu K \ 

BF FfT o:S H “ §7/ 

(bv o „ 

bS a 1 _ / bw H bo„\ {bu , bw,,^ 

bl yR ^^Kby SF)~ u °^ a - ( 8i~ i W) 
Therefore, the entropy field can be expressed as 


(3) 


bS _ bs 0 * n bS,-) 
bx 'BF~ + ' 2 r' * 5a: 


bSj, 

by by " by 


y 


bS 

bz 



bS a 

bz 


J 


(4) 


Assume that each coefficient a t ... a* is constant in the 
entire region of the flow field where it is not zero and is small, 
so that terms of the order of a, 1 or higher can be neglected. 
In this flow a basic flow field exists, represented by the velocity 
vector V 0 (wo,Cb) and by the entropy distribution S 0 , on 
which a linearized flow is superposed, represented by a 
summation of N three-dimensional flow fields, each of which 
is proportional to the corresponding coefficient a n . The 
basic flow field is a two-dimensional flow if Cartesian coordi- 
nates are used or an axially symmetric flow if cylindrical 
coordinates are used and can be determined by known 
methods. In a similar way, a general three-dimensional 
flow field can be assumed for basic flow if the flow field cap 
be obtained by simple analysis. 

The equation of motion obtained from continuity, momen- 
tum, and energy equations can be expressed in Cartesian 
coordinates in the form 

bu{. _t^\,bp/ t£\ , bw/. w 2 \ uv ( bu biA 

bx\ a 2 / by \ aV + b2\ aV a 1 Vby^ bx) 

uw{bu . bw\ wv {bv . bw\ . ... 


while, in cylindrical coordinates, the equation becomes 


For constant stagnation enthalpy in the flow field, the rela- 
tion between entropy and rotation states 

curl V X V grad S (2) 
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Bv use of equations (1), expression (5) becomes 
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-M 1 - 

Uo J \ , 5 po/, 

'atrby V 

Po*\ 

■“s? y 

u 0 Po/bUo , dgo\~j 

no } \&y dx/J 

U Q = Fq cos <?q| 

O) 


, du, (. a/ 
*" bx \ a?. 


\ Pfl*\ UgVofbo, , bu,\ 

„ a 0 v Uo s \bx + by) 

Po= Vo sin tpo ) 

and, because higher-order terms are neglected, 


bug (' 2u a u M +u 0 1 A n '^ boo(' 2ooP 1t +Oo t A K '^ _ 


dx \ 


a* 


dy V 


at 


> 


(dug ago \ ifog 0 xl,+itor»+gotfn ~[ 

\dy + dx/ do 1 J 


( 7 ) 


where 


y ^ 

-4,=-^— ; r (now^+ror.) 
Go 


The basic flow is a two-dimensional flow, and the first part 
of equation (7) must be zero; therefore, the second part of 
the equation must also be equal to zero and, for each value 
of n, 




bz 


du*/. «A, »A 'MqPq /do» | au A 

dx V by \ a 0 V a* V dx ‘ by ) 

_da o / 2 'u a u n + u a 1 A x \ , Sr o / 2p 0 g,+g 0 M, \ . 
dx V. a Q 2 y by \ «o l / 

/du 0 . do 0 \ 

\djA dx,J 


« or 0 xl,+UoP«+Potf, 

ffo 1 


( 8 ) 


Equation (8) is a differential equation of hyperbolic type 
because aa*+Po*>Uo* and the characteristic surfaces are 
cylindrical surfaces perpendicular to the plane z— Constant 
with generatrices coincident with the characteristic lines of 
the two- dim ensional flow. This can be seen from the fact 

that the coefficients of the derivatives and 

of equation (8) are expressed as functions of the properties 
of the basic flow and are the same as the coefficients of the 

bo o 


■t which define the basic Sow, 


, . ,. du 0 bo 0 j 

derivatives -5- — * -5 - — 7 and t — 1-3: — 
ox oy ay ox 

and the coefficient of is one (see, for example, reference 4, 

page 282). Therefore, the disturbance flow field V H (u K ,r x ,w K ) 
can be obtained by the method of characteristics by moving 
along the characteristics surfaces which are cylindrical. The 
characteristic net, which for the general case must be drawn 
in spatial coordinates, can be drawn in this case only once 
for any value of a and n and is equal to the net of the two- 
dimensional flow. Equation (8) can be transformed for 

practical use. In the equation, the terms in UgVo, -~^J ! 

and ^ are known terms and are given from the two- 
dx by 

dimensional flow field. 


Assume, in the plane xy, polar components for F and <p for 
the velocity, defined by 


where 


a*= F* cos <£o — <p*Vq sin <pA 
r»=F, sin ■aj+pJ’o cos <po ) 

F=F 0 +Sc.r,I 

p=?i)+S ll »?i / 


( 10 ) 


(11) 


If n is the normal to the projection of the streamline in the 
plane z= Constant, the derivative ^ is expressed by 


bS 


bS p 


bS u_ 
"dy V 


bn bx V 

By use of equations (4) and (1), equation (12) becomes 


Cl 2) 


bS = _f' 

bn V dx 

'dS, 




\ Pq+ UmP it I 

/ Fo+Sa.FA 


/bS a , dS*\ ao+ 2>*u» 
U*/ +2Lj m ayAro+Sa.F; 


(13) 


By considering only the Iowest-order t erm s, equation (13) 
becomes 


on r orq 


bS> 


s«.[- 


bng 


-^(p»F,-p 0 F,)+ 


(14) 


But, 


- (r* Vo- Vo F,) + ^ («» Vo- Uo F„) 


by 


^■tfbSo u 0 , bS„ o 0 \ n 

=-<?Ao (j^v+ W tJ =0 


because the right-hand term contains as a factor the variation 
of entropy along the streamline of the basic flow; therefore, 


dr* 

bu m _ 

Chs* 

bS* / 

dx 

~ *>y ’ 

-yRVo 

da 0 V 


<-> 


where n 0 is the normal to the streamline of the basic flow in the 
plane xy. 

After several transformations the following equations can 
be obtained (see reference 2): In the plane c=Constant along 
the characteristic line defined by 


^|=Xi=tan (jSo+^a) 


(16) 
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the following equation is valid: 
dw n 1 sin jS 0 ten p 0 . 1 dV % 


dz F 0 cos (/3 0 +?o) F 0 dx 

sin 3 (So 1 dS, 


-ten f ? 0 

F, 


cos fa+/3 0 ) yR dn 0 +<P ” B ^ + \* a ° l ~ 0 (1 7) 

and along the characteristic line of the second family defined 
by 


^=X 3 =ten fa-j8c) 


(18) 


the following equation is valid: 
dw n 1 sin iS 0 ten 0 O . d F» 1 


dz Fo cos fa— /3 0 ) 
sin 3 /S 0 


dx Fo 1 
1 dS ; 


tan fi 0 d -^~ 

—‘pnB i -{-~- C%=Q 

Vq 


(19) 


cos fa— (3 0 ) yR dn o 

where Bi, B a , £?,, and 6V are coefficient functions of x, y, F 0 , 

<p 0 , and a ^ on & Ibe characteristic lines of the 

first and second families at each point and are independent of 
F„, <p n , and w n , and, therefore, can be determined once for the 
basic flow and used for any kind of disturbance flow in the 
limits of the approximation accepted. The coefficients B u 
B t , Cu and C» are given by the following expressions: 


Br- 


1 


COS do cos fa+ do) 


B 2 = 


COS do COS fa 


T sin 4 /3 0 

dS 0 

±( 

dV A 

cosfa— /3o)"j 

7)1 yR 

dn 0 

FoV 

. dx /: 

, sin /S 0 J 





(20a) 

f sin 4 /3 

o dSg 

i 

( d Fp 

\ cos fa+doTj 

L yR 

cItiq 

Fo 

\ dx , 

A sin /3 0 J 





(20b) 


°‘-r, ;(w),fa‘ , « + 

j~ 2 dS 0 sin 8 do 
sin 3 do 


7—1 


2 sin 3 do cos 2 do) 

1 / d FA cos(y 0 — do) ~ 1 

7 R dn a cos(y>o+do) Fo v dx A cos fa+ 0 O ) cos 3 doj 

) . ( 200 ) 




7—1 


2 sin 2 do cos 2 do. 




2 dSg sin 3 do . 1 / d V a \ cos fa + do) 

7 It dna cos (tpo do)” F 0 \ dx A cos fa— do) cos 


[ 

( 1 + 2sIn 2 do) 


J’do] 


(20d) 


where is the derivative along the characteristic of the 

first family and (^jr) , along the characteristic of the second 
family. 

In order to determine the value of w K at each point of the 
characteristic net, the following relations can be used in the 


plane 2 = Constant: 


_bw 

bp 

~*>y~ 

dz 

_bu 

bw 

~bz 

bx 

_bv 

bu 


( 21 ) 


dx by J 

and, if s 0 is a streamline projection in the plane 2 =Constant, 


bw bw u , bw v 


uiv irw u, ,uuj v U . U , u » 

V + ~by V V v + V* + bz 
In the approximation accepted, 

3F „ bV> 
bz ~^ a '~bz I 


bV bV ,bS a* 


’ &2 bz ylt\ 


» ( 22 ) 


therefore, 


bS ^ bS n 

bJ-So.-sj 


1 bw n 1 5F» , bS n ti 0 3 
Fo bs 0 F 0 bz + bz yR F 0 * 


(23) 


(24) 


where s 0 is the streamline of the basic flow. 

Equation (24) permits the determination of the value of 
w at each point of the plane 2 =Constant as a function of 
the local variation of F„ in the 2 -direction. Therefore, by 
means of a step-by-step procedure, all the flow field can be 
determined by working in planes z— Constant where only 
one characteristic net must be used, by means of equations 

(17) and (19), and by determining the value of ^ at each 

point of the net by means of equation (24). The calcula- 
tions are started from the flow field defined, or along a sur- 
face which is not a characteristic surface, or along a charac- 
teristic surface and a stream surface. The flow field at the 
starting surface must be determined from the boundary 
conditions. If the starting surface is a shock wave, the 
flow at the shock surface must be obtained from the physical 
properties of the shock wave related to the boundary con- 
ditions considered. Relations between boundary conditions 
and shock waves are presented in detail subsequently for the 
problems considered. 

Similar equations can be obtained by using cylindrical 
coordinates in place of Cartesian coordinates. 

Equations (17) and (19) remain the same. Only the first 

term changes: — becomes and the coefficients B, and 


B a become 

Bi = 


1 sin ftp sin fa— /S q) fa 

cos £ 0 cos fa+£o) V 1 


□ /__ 1 sin ffpsin fa-H?o) . „ ( 

2 cos jS 0 cosfa— /So) V 1 U3 


( 25 ) 
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Equation (24) for cylindrical coordinates becomes (see 
reference 2, equation (8c)) 

/£> _ ar, , a 0 3 dff, to, sin 8 a . . 

V }t — Contfartf y 50'" yRV a y d0 y ^ 

The use of the linearized characteristic system is simpler 
in many cases than the complete characteristic system and 
reduces in some practical applications the extent of the 
numerical work required, especially if solutions of many 
similar problems are required. The same concept can also 
be easily applied in the field of nonsteady or relative motion 
of flows. 

The flow field around slender bodies without axial sym- 
metry can be obtained by means of linearized methods, and 
some practical three-dimensional problems not analyzed 
before in the approximation of nonlinearized flow can be 
analyzed by this method. 

In the next sections some typical possible applications are 
presented. 

SOME TYPICAL APPLICATIONS OF THE LINEARIZED 
CHARACTERISTICS METHOD 

TWO-DIMENSIONAL ROTATIONAL FLOW FIELDS 

Two-dimensional potential flow permits hedograph solu- 
tions, and, therefore, any kind of two-dimensional super- 
sonic potential-flow solution can be obtained in the hodograph 
plane; the numerical solution in the physical plane then 
requires only the construction of a characteristic net in 
order to find the position in the physical plane of each 
point of the hodograph plane. The solutions of problems 
in which boundary conditions are given only along a stream- 
line are very simple in the approximation of potential flow 
because in this case the velocity is constant along a family 
of characteristic lines (single-wave flow); however, similar 
calculations for rotational flow are much more involved, 
because a step-by-step procedure is required for the solution 
on the hodograph plane as well as for the construction of 
the characteristic net. By means of the linearized charac- 
teristics method a rotational flow field can be considered as a 
modification of a potential flow field, and the linearized 
superposed flow is the flow which takes into account the 
effect of the presence of shock waves and the effect of rotation 
in the flow. 

Consider, for example, a two-dimensional profile which 
produces a shock wave at the leading edge (fig. 1). If the 
profile is curved, the shock is. curved, and the flow behind 
the shock is rotational. Assume that the flow field behind the 
shock can be expressed as 

«=u Q -fu I 

r=r, +r, 
s=s, 

where u n and v a are the velocity components defined by a 
potential, flow field, which in this case is a single-wave flow, 


C=E> a + C 1 'j - 
<P=<Po+<Pi)' (27) 
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Fic.I'ee 1. — Application of the linearized characteristic system to two-dimensional rotational 

flow. 



the characteristic lines of the first family is 


zero, and u lr c u and S L represent the flow field that takes into 
account the reflections occurring at the shock and the effects 
of the entropy gradient. The flow field represented by V a 
and <p a can he immediately determined, and T u and <p 0 are 
constant along the characteristic lines of the first family, 
which are straight lines. The characteristic net of the flow 
V a>Vo can be drawn in a short time. Then, along the charac- 
teristic of the first family, 


dV i -l- tan Ci dx^O 


where . 
and 



1 1 1 

fdV o> 

1 cos {vo— 0o) 

cos 

00 COS (p o +0o) To ' 

i dx j 

sip 0o 

/dVo' 

\ 1 COsOpo— 0o) 

1 

(l 1 7-1 ^1 

\dx , 

'i To cos Opq + 0o) cos 2 0o 

'n 2 sin 2 0o/ 


(28) 

(29) 


and, along the characteristic of the second family, 
pr dFi+tan @o d<pi-{- dSi~b 

0o+^ — ^-rw— — ry) = 0 (31) 
V a I o \ 2 sin 2 j3 0 cos 2 0 O / 

Equations (28) 'and (31) can be simplified by' introducing 
the gradiants (Vi)* and («j+0o)* along the x-axis defined. as 




(32) 


It can be seen from figure 1 that, at any' point A of ordinate 
y A on the characteristic line a, crossing the axis at A 0 , the 
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the form 


variation along the characteristic line of the second 

family at A is 

r/^l_ fn r cot (go+foW l _ .1 

sin 2 (^o+ft).i 0 
(33) 

because ft and <p 0 are the same at A and at A 0 . Therefore, 


WdVo\ = 


a 


Vo \dxJ A 1 — by A 


(34) 


where a and b are constants along each characteristic line 
of the first family. Then equation (28) along the charac- 
teristic line of the first family becomes 

dVi , „ „ , sin 2 ft JO 

-p-- — tan ft d<pi-\ — pg—Goi— 


cos(po— ft) 1 ady 

lpl sin (vo+ft) cos ft sin ft 1 — by 


+ 


, 

Kv 2 sin 2 ft 

Therefore, if 


7—1 \ a dy cos (lo—ft) 1 

J 1- 


y • ' ■ -j 

T>y sin (<o 0 +ft) cos* ft 


=0 (35) 


.2 (PA 1 

“ sin 1 (#><,+ ft) Vo 

i _ (ffo+fth 
sin 2 (e’o+ftj 


(36) 


the equation along the first characteristic line becomes 
tan pa dSi + ^i -j log (1 by)— 

( 1+ i (37> 

and the equation along the second characteristic line becomes 
d Fl tan fig d Vi + - ^ dSi+ ' 


V„ 


Vo Vo ( Un2 Po+ 2 sin*ft cos 1 ft) 0 (38) 

All the copfficients of equations (37) and (38) are constant 
along characteristic lines of the first family and can be 
calculated at few points on the z-axis. 

The coefficients of equations (37) and (38) are independent 
of <pi and Vi; therefore, with one calculation from points A 
and G of the net of figure 1 the values of Vi and <p t at a 
point B can be obtained directly without the necessity of 
an iteration process, with terms of the order of (Az) 1 also 
included in each step. Indeed, if all the quantities which 
are variable along the characteristic lines are expressed in 


( 39 ) 


then 


■'=«'+(sl 1 '+§T+« 

(siHsCl+S^+W 

t40) 

Consider now equations (37) and (38). At point C equation 
(37) has the form 

(^X-(^)o‘‘' + (w)c 6c+, '«(3f)c + P, «(e)c 

while at point B it lias the form 
then, from equation (40), 

2(V,-F, c )=(^)y+(^)y 

X+ P -(s)> < 4i > 


'd( 

M asr 


+<pi 


'd<_ 
*\dx) B 


+ Vl c(% 


But, 

Indeed, 

and 

[( wX ^+ Cfe 1 )/*] 41 

and also 


Ax-i- 




Therefore, equation (37) along the first characteristic fine, 
terms of the order of (x B —x c )* being neglected, can be written 
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in the form 


0 V~ 1*) r\ { vib ~ 'P'c) tatl &*+ 

* 0 r 

^ sin* *,+| (log '&*) [(*.+,■*)- 


tan 0„ - ^1 • 


7 — 1 


2 sin* /So^/J \6/c 


?) =0 


(42) 


and equation (38) along the second characteristic line, term 9 
of the order of (a : s — 3 being neglected, can be written in 
the form 

0 V~ *\t) ( ¥ ’ i «~ ?, u) ( tan ^o^+tan /J,J + 

S \fT A ^03+ sin-’ /3 o 4 ) + 

f [t? ( tan * p0+ 2 sin*0 o cos* (9 0 )1 + 

[fo ( tan ‘ J3o+ 2 sin*0 o cos* J]J ( 1 ^ ^) = 0 

(43) 


procedure (for example, reference 3), and in a region inside 
the body as shown in figure 2 in order to determine all the 
fiow field necessary. Thus, a reduced number of points of 
the characteristic net is chosen at which the superposed flow 
field for each boundary condition different from the basic 
shape will be determined. The number of points required 
depends on the magnitude of the superposed flow field; how- 
ever, the 1 number is usually small, because for each step 
between two points A and B the disturbance velocity can be 
expressed in the form 

(F0-(F0.+(^)/x+(g-')^. (46, 


where the term 



is also included because, in the 


differential equations (44) and (45), the coefficients of the 
differential equations are independent of the solution and 
are known at both points. If the entropy term Si is neg- 
lected, by applying finite-difference methods the velocity 
components at a given point g5 of the characteristic net can 
be expressed from the values at two points f5 and g4 in the 
forms 


Because of the possibility of considering directly in the 
calculations terms of the order of (Ar)* also, large steps can 
be used in the characteristic net and the effect of entropy 
gradients can be easily evaluated after the basic character- 
istic net and a few streamlines of the basic flow have been 
determined and the coefficients of equations (39) and (40) 
calculated at a few points on the axis. 

All the coefficients are constant along each characteristic 
line of the first family; therefore, in going from B to D, only 
the terms containing Si and y must be changed and the 
calculations are simplified to some extent with respect to the 
rotational-flow characteristic calculations. 


AXIALLY SYMMETBIC PLOW FIELDS 


For axially symmetric flow, the equation of linearized 
characteristics becomes 


1 dF, 


-tan 


dp, .sin* 00 dSi 


r. 


yiT d^ l5/+ Fo Cl - 0 (44) 


To dx R H0 dx 
along the characteristic line XL=tan(j3o— p 0 ) and 

l dr, 


To dx 


*“ ^ +si y 7r-**'+ v. < 45 > 


along the line \ 2 =tan(p 0 — 0o) where Bi, B/, C\, and Cj are 
defined by equations (25) and (20). 

The introduction of equations (44) and (45) simplifies 
noticeably the numerical calculations without affecting sensi- 
bly the precision of the results. The practical use can be as 
follows: A basic body shape is determined first by means of 
characteristic calculations, and the characteristic net is then 
obtained. The basic calculations must be extended in a 
region in front of the shock wave, determined with the usual 


(TU s =(l’,) /s J~+(p I ) /s f±5 +(¥Ji)ts i_P. (48) 


tei)«= 


1 


m— . g— p 
d+n + 1 + r 


>*>4 
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+ « 


l)ra T+>"" {p,) * t T+7T ' 


(49) 


where (p,) ft , (<pi)rs, (T',) t i, and (Idj/s are functions of the 
boundary conditions considered, while the coefficients L, 
m, n, p, q, and r are functions only of the basic flow field 
and, therefore, must be determined only once for any 
boundary condition considered. These coefficients are 


jr = (Fo tan- flo)«*+( T r a, tan ff„) fS 


«“(**-* J K VoUB^+iVcUBSU | 


n=(x l5 - J : ft )[(C) I 4+(C,) ts l j 

__ ( T r o tan ggh.-KV, tan 0ohs 
V 2 


(50) 


3 — x /s) K I r o)*e(fi/)*s+( T o )/s(5i0/*l 'g 


r=(-c fS -^)[(C' i ) l5 +(C' 2 ) /t li 


If the entropy terms in 


dSi 
dn o 


are considered, two more terms 



940 


REPORT 1102 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


in (Si) t6 and (Si ) f4 must be considered in equations (47) and 

(48) and two more terms in (S^/t and (Si) ti , in equation 

(49) . However, because these terms are small, they can 
usually be neglected in practical calculations. 

Because all the coefficients are determined only once, the 
determination of (Fi ) t5 or for each boundary condition 
is simple. Points at the boundary can be investigated by 
means of equations (47), (48), and (49). Points on the shock 
can be analyzed in a similar manner. A practical calculation 
can be performed in the following way: The basic flow field 
and the number of points in which the superposed flow field 
will be considered having been determined, the new boundry 
conditions (shape of the body) are placed in the character- 
istic net (fig. 3). If A is the point where the basic body 
departs from conical shape and OBC is the new boundary- 
condition, the flow field between the surface of the conical 
body OB and the conical shock OE are known from cone 
calculations; therefore, the values of F and tp at each point 
Ba, la, 2a, 3a, and Ea of the basic characteristic net are 
determined, and by difference the values of Fi and <P\ can be 
obtained. From al the flow at the point L of the boundary 
can be determined from equation (48) where yn is known at 
Tv, and the coefficients at L can be determined by- linear 
interpolation between hi and cl. From the values of F, and 
<pi at al and Ll, the corresponding values at bl arc inter- 
polated. Then all the values for the line b can be obtained. 
For the determination of F, and y>, at a point F on the shock 
wave, the equations of the shock wave and equation (47) 




FIO0RE3 —Practical application of the linearized charaeteristlcssyateTn to axially symmetric 

flows. 


are valid. At the point F the values of and F a arc known, 
and from the equations of the shock wave the value of 


as function of y> can be determined. (The value of y> fixes 
the deviation across the shock wave.) Therefore, 


(F)r=(F 0 )r+(Fi)r=(Fo)r-f 



(^i)r 


(51) 


Then, equation (47) applied between the points 4b utul F 
gives 

Wr=(lPl Mr' 

*o 


V 

=(F,)* 


1 — n 
1 


'(<Pl)n 


m-\-L 


t+7 


m—L 


(52) 


and the value of (<pi) r can be determined. 

The work required in the calculation of the flow field for the 
basic body and the determination of the coefficients L, in, n, 
p, q , and r can be reduced to a minimum if conical bodies are 
assumed as basic bodies for the calculations, because in this 
case the basic flow is available in tabulated values (reference 
5) and the coefficients L, m, n, p, q, and r arc functions only 
of the polar coordinate ifx. 

For conical flow, the calculations can be performed in the 
following way-: From conical-flow calculations, the values of 
D r and i\ as functions of \p are known, where o r is the radial 
component and v K , the normal component of the flow field 
in polar coordinates referred to the limiting velocity (fig. 4). 
The following expressions for F 0 , /9 0) and <pn can be determined 
from the cone calculations: 


F*= t>* 2 + D r 2 

1 i 2 . 

vo—i+y 


tan 7 = 


(53) 


Therefore, 0 O + <Pa and y> a — /S 0 are known as functions of 
From a point A on the conical body, the characteristic line 



x 

Fioure 4. — The basic net when conical flow k assumed as tba basic flow tor the linearized 
characteristics method. 
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of the first family AC and of the second family AD can be 
* drawn. The line is defined by the expressions 

Vs— yx ^ tanUg+^+tanCff-f 

£k—x.< 2 

Vr- y.< _^ tan(/3— p^+tanOS— <p) r 

x A —x F 2 

In order to construct a characteristic net that requires only a 
small amount of calculation, the points of the net are chosen 
along straight lines from O so that the conical property of 
the flow can be utilized. The net can be constructed by 
fixing the steps along the body. 

When point B is chosen along AO, the point E along AC 
is determined by drawing BE parallel to AF. From E and 
F (along OE and AD) the lines GE parallel to LF and FG 
parallel to EH can be determined, and the point G can be 
obtained. From G and H, point M can be obtained, and by 
proceeding in a similar way, all the characteristic net can lie 
determined. The coefficients of equations (50) are the same 
for each point N, F, and E along the same radius. The 
tabulated values are given only in the region between the 
body and the shock wave; however, the calculations can be 
extended by means of conical calculations. For example, 
from the following equations (see, for example, reference 4. 
p. 243, and the following pages): 

(/lu.-r — 


Cr+C* COt \p 

_ 1 (y — i k i ~ ^ vj _L 


(C.)*,. +4 *=<>«)*£. cos A^+(jB— c r ) t /sin A^ ) 

/ (5 5) 

(prV c+ A#=(0* c sin Atf— (fl— v r ) fc cos A^-f(2?) #c ) 

(where R is the radius of curvature of the streamline in the 
hodograph plane); therefore, the characteristic net can be 
extended to the outside flow. For conical flow the coeffi- 
cients B\, Bi, Ci, and C% can be determined as functions of 
ip. For conical flow, coefficient B/ of equation (25) becomes 


l sin go sir 

cos 0 O cos (<fo+Po) L ; 

1 d Vo / djV\ cos (yo— j3 0 )~ | 

T’o dip \dx /x* sin 0o J 


sin /3 0 sin (<po— 0 O ) 


sin /So 


and Ci becomes 


. , _ 1 d To / 

1 To dip War 


(fifiXt cos'/Sq (® ia ^ 0+ 2 7 sin ! /S 0 ) 


1 dV*/dt\ cos (<po — /So) 1 A , T-l \ 
To dp \dx)\ 2 cos fo)+j3 0 ) cos*j3 0 V ~ r 2 sin 3 j3 0 / 


rl£) = 

v dx A, 


sin ip sin (<?q + /3 0 — p) 
V cos(<po+/So) 


/ dp\ _ sin p sin(i/' — pqT/3q) 
Wx/xj - y cos(^o— /S 0 ) 


u=r Bi r 0 (-jr) ] 

L vtyAiJ* 

1 |~ r sin/Sp sin(y 0 — j3 0 ) , 

cos 0o L V ° sin p sin(<ft + /3 0 — py 

fdVo\ 1 sin (i/>— yo-r/S 0 )~[ 

\ dp )+ sin 0 a sin (<po+/3 o—PU 

K* r *(w\l 

_ 1 r y sin /So sin (yo+/3 0 ) 
cos/SoL 0 sin p sin {p— pa+jSo; 

/ d V 0 \ 1 sin (<po+<3q— j') '! 

\ dip ) t sin 0o sin (p— po+/3o)J 

= To (jdf\ [co?fr ( Sm ^ 0+ 2W3“ 

1 /. | T — 1 \ sin (p— <pg+ ffo)~ j 

cos 3 0 o\ ' 2 sin 3 00 ) sin {<Po+ 0 o—Pu 

~T,(cdf\ [coiFS ( sm ’ |S,+ 2"am* (ij- 

1 , y — 1 \ sin(yo+/3 a — p) ~\ 

cos 3 0 a \ ‘ 2 sin 3 0 a ) sm(p— vo+/8ii)J 


r r +r, cot ip 
2d * 


dip Vo T 2p»* (58e) 

(r-l)(l-To 3 ) 

If equations (58) are used, the coefficients of equations (50) 
can be expressed as follows: 

1 'i 

™=4 OAs— ^,rt)(a« 4 +a*s) 

71 (Pis-Ptl) (c It + c, s ) 

i r (59) 

g=4 (Pu— \As)(As+&**) 

r=~ ('P t f,— pfs)(d !S +d gS ) 

* 

Because the coefficients a, 6, r, and d are functions only 
of P and of the free-stream Mach number M, they need to be 
calculated only once for different values of M and given in 
tabulated form; therefore, the calculation of any flow field 
for which the basic flow can be considered a conical flow can 
be reduced to the solution of a few linear equations with 
known coefficients. 



942 


REPOHT 1102 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


In order to try the method, the flow around an ogive, as 
shown in figure 5, has been determined at M=3.016 by the 
method of characteristics and by the method of linearized 
characteristics. In figure 5 the usual characteristic net is 
shown, while in figure 6 the linearized characteristic net and 
the basic body are presented. As a first basic body the 
cone chosen is the cone tangent to the apex having a cone 
angle of 12.5°. The cone chosen is not the most convenient 
because the values of <pi at the end of the ogive are large, 
and a better approximation would be obtained if a cone of 
smaller cone angle would be considered as the basic body. A 
cone somewhat different from the ogive considered has 
been chosen in order to have some information on the 
approximation of the method for sensible variations of the 
shape of the body from the basic body. 

Downstream of the characteristic line AB, the velocity 
components F, and <pi would become large because at the 
surface of the body (lie component <p is quite different from 
the component po of the basic flow; therefore, the flow deter- 
mination in the region downstream of the line AB has been 
considered as a new problem, defined by the flow along AB 
and from the streamline that represents the body shape. In 


this region a new basic flow has been considered. Again, a 
conical flow field has been assumed as the basic flow. From , 
the values of <p and V at B and A, and from the order of 
magnitude of <p downstream along the body, a conical flow 
field that would give a small disturbance component in this 
region has been selected. The cone chosen for the second 
part is a 5° cone at M=3. 077. The cone is entirely contained 
within the body considered in the region used in the calcula- 
tions as shown in figure 6, and the most convenient region of 
the conical flow field is used for the calculations. 

In order to pass from one basic body shape to the other, 
the. components F, and <pi along the characteristic line AB for 
the second basic body must be determined. This operatiQn 
can be performed by determining the value of TV-f-TV and 
at the points a, b, and n of figure 6 for the first basic 
bod}', by interpolating the values at the points a', b', and n' 
between characteristics 1 and 2 along the characteristic of 
the other family, and then by determining the new value of 
TV' and <pi" for the second basic body along AB from the 
expressions 
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and 

r=TV+TV 

'?=<£<> -tVi' 

where Vo" and ^o" are the values of V and <p for the second 
basic body at the position considered and V 0 ' and #/, the 
values for the first basic body. 

The calculations by linearized characteristics required the 
solution of 11 linear equations of the type of equation (47), 
(48), or (49) with one unknown and the interpolation of four 
points, which can be done in a very short time (of the order 
of 1 hour) when the net is drawn and the coefficients of the 
equations are determined. The pressure distribution ob- 
tained is presented in figure 7 and is compared with the pres- 
sure distribution obtained from the exact method. 

For the back part of the body where the flow differs slightly 
from parallel flow, a cylindrical body with uniform flow at 
different Mach numbers can be considered as the basic bod}'. 
In each region of the flow the Mach number for the basic 
flow is constant; however, the component u can be approxi- 
mated conveniently by changing the basic-flow Mach 
number. 

When the body has a tail, a conical solution as proposed in 
reference 6 for flow inpide a tube can be assumed as the basic 
body. In this case, each streamline of the conical solution 
ean be considered as the shape of the basic body (fig. 8), and 
the flow field can be obtained from conical-flow calculations 
which can be determined easily in the hodograph plane. By 
changing the strength of the final shock of the conical solu- 



Fmvu T. — Comparison between results of the calculation by the characteristics method 
and by the linearized characteristics method obtained at .\f-8.018. 



Figcix 8.— Conical flow for tail of axially symmetric body. Initial .Yf-3.2. 


tion, different ratios between maximum cross-sectional area 
and tail area can be obtained. Any part of the streamline 
can be assumed as the basic body shape. 

The present method does not require the existence of a 
linearized solution and, therefore, can be applied also at 
bigb Mach numbers. This method permits obtaining the 
shape of the shock wave and taking into account entropy 
variations. High precision can be obtained by using several 
basic flow fields for the different regions of the body con- 
sidered. Because of the simplicity of the calculations, the 
systematic - calculations and tabulation of coefficients of 
equations (58) for different cones and different Mach num- 
bers would be of great practical interest. 

Tabulated values can be obtained also in the following 
way: For each cone OC of cone angle ^ considered, a super- 
posed flow must be calculated as shown, for example, in 
figure 9(a). The values of V l and (?i ■ for this superposed 
flow are obtained at given points of the characteristic net. 
Because of the linearization of equations, if the superposed 
flow changes in intensity, the values of Ty and at every 
point change proportionately. Because of the conical prop- 
erty, if the point A (fig. 9 (b)) moves along OC, the flow 
field changes in scale; therefore, the effect on a point E due 
to the superposed flow field Ap t starting at A is equal to the 
effect of a linearized flow field starting at A' and of intensity 
A«?i at a corresponding point E' defined by 

X*'... Xa' _V3’ 

x a Xj. Vr 

Therefore, when the flow field for the disturbance AB is 
determined, the effect of any disturbance of the type of the 
disturbance AB in the entire flow field can be obtained. 
Then any body shape can be considered as a superposition 
of flow fields of the type of flow corresponding to the dis- 
turbance A <p placed along the basic cone. From the simple 
calculation of the flow for the shape AB, the velocity can be 





B 


i 

8' 

t 

\ 

A' 

(a) 



(a) Disturbance distribution. 

(b) Position of the disturbance. 

Figcm 9.— Determination of linearized flow of constant Intensity superposed on conic a 

flow. 
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determined by means of the equation 
V=V 0 +±V n 

i 

where n is the number of the superposed flows that affect 
the point considered. 

Consider, for example, figure 10. Several linearized flow 
disturbances must be superposed on the basic conical flow 
field. First, a superposed conical flow at 0 that can be 
obtained from conical-flow calculations must be considered. 
At Ao a superposed disturbance must be added in order to 
satisfy the boundary conditions at A. If the calculations 
have been performed for the disturbance at a^ and for 
£^=1, the velocity at A can be obtained from 

F^fF^ + fFOx+GYx 

where (F 0 ) A is the velocity of the basic flow at A, (Fj) x is 
the velocity disturbance at A due to the conical flow super- 
posed at O, and (F*)a is proportional to 


<PJ a=VA~ [(%) A + (?>i) a] 


and can be obtained from 




where b is determined by 


Xag Xt y„ 


In a similar way F 2 is determined at the points B, C, and 
■so forth. At B another linearized flow field having velocity 
components V 3 and <pz must bo considered where 


( <Pi) b — <Pb — [ ( IPo) b + ( (Pi) b + (ps) b] 

and 

(F S ) S =(F,)>^ 

m)o 

where b is defined by 


x^ x» y b 



Viouri 10.— Application of the linear lied characteristics method having constant 
disturbance superposed on & conical basic flow. 


Because of the rapidity of calculation, the variation of any 
geometrical parameter can be investigated in practical appli- 
cations without the necessity of a large amount of numerical 
work. 

CONICAL FLOW FIELD WITHOUT AXIAL S1TMMETRV 

The calculation of slender bodies without axial symmetry 
requires the determination of conical flow without axial 
symmetry, which can be done by means of the linearized 
characteristic method. The basic problem of the de- 
termination of conical flow consists in determining the shape 
of the conical shock wave produced by the body. Wien the 
Bhape of the shock is determined, the flow field around the 
body ctln be obtained by means of numerical calculations 
(see, for example, reference 4). Because the relation be- 
tween the shape of the body and the shape of the shock wave 
is not known a priori, the method of linearized characteristics 
can be particularly useful for flow' determination of this kind. 
An approximate shape of the shock wave is assumed as the 
basic solution and the flow field inside the shock is de- 
termined; then a linearized flow field is superposed in order to 
satisfy the boundary conditions at the body. The calcula- 
tions are simple if the basic flow can be determined analyt- 
ically or numerically without a large amount of calculation. 
For example, for slender bodies the basic flow can be the 
axially symmetric flow' for which values arc available in 
tabulated form. Consider a conical shock wave which can 
be defined in polar coordinates as 


cos nfl+S(&»), sin mO (GO) 
a l 


where all the values of \[> n and are small so that terms of 
the order of \j/ n * can be neglected. Such a shock wave is 
approximately of circular cross section, as is found for slender 
conical bodies. If the flow is assumed to havo a symmetry 
plane, the second summation of equation (60) is equal to zero. 

The velocity components in the radial direction (e>)i, in 
the tangential direction (e r )i, and in the direction normal to 
the shock (tv) i in front of the shock wave are (see fig. 11) 

(Pr)l = Fi COS \ 

Mi= — Vi sin cos cc > (61) 

(Pr)i= — Vi sin \p,sina) 


W'here all the velocity components are referred to the limiting 
velocity and a is the angle of the dihedral between the plane 
normal to the shock wave and the plane containing the ref- 
erence axis. Across the shock wave the following relations 
are valid: 

(iv) I C4?v)s=^Y[l-(» r )i s -(t? r )[ i ] 1 


^ p~ -(l — V\ cos 2 F ( 8 sin 2 sin 8 a) 


{'V n)i — 


7—1 (1 — F t 2 cos* jy— T V sin 8 sin* 
7+1 Fi sin \p, cos a 


«) 


r(62) 


(Pr)i — (Pr) 2 
(t>r)l = (t>r)2 
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If the velocity components behind the shock v Tl and w of 
the flow field in polar coordinates are considered at each 
point &=&),+ cos n 8 (fig. 11 ), 


(r ri j^=tv= cos p t 


(63a) 


. . 7 — I (1 — Ft* cos 1 p , — Ft* sin 2 p, sin 4 a) 

7+1 Ft sin p. 

Ft sin p, sin* a (63b) 


(w)+ t — —F t sin p, sin a cos a+ 


but 


7 1 rr J 4 , • • 4 i » ] \ SID d 

— r-T (1 — V 1 COS* p, — 1~ sm* p, SUL* a) ff — « ; 

7+1 v r r ' Ft sm p, cos a 

(63c) 

dp np x sin.nfl 

tan a=a= . * = >) — — +— 

sm p dg sm p 


Therefore, the velocity components behind the conical shock 
wave of equation (60), if terms of the order of p x are neg- 
lected, are 


(Cr)*,=F cos (p Q ),— Vi sin (*0.S(*.). cos nd 

\ _ r~ 1 1 — ly-cos* (pg), 

K) *' 7+1 F sin (^, 




7 _1 
7 + 1 


cos n ' 



7—1 

7+1 


/ l— F* cos* 

V Ft sin* ^o, /A 2 - 1 


n(P»)t sin n6 



Then, the velocity components behind the shock defined by 
equation (60) can be expressed as: 


f>r)# = (t>r „)*+ X COS n6\ 

( v „) # = (+,,,) i +^(p x ).(v ,)« cos nB\ (65) 

(u% = 2XW. w.n sin nd ) 

where iv 0 , (rj », (cj*, and ic, are independent of d, p, a and 
l\ correspond to the flow field for circular conical shock 
Pa,, and (r,)„, (r,)„ and (w)* at the circular cone p=p 0 , are 
independent of n and are defined by 


[(c r )»l* 0 ,= -Fi sm K-O«o)* 0 , 

w-k-vMJ&X 




and 


( 66 ) 


/ Sp «o\ 

— 

*>■„+» »o cot p 0 - 

v * 


a 0 3 
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It can be shown that the flow field defined as in equations 
(65) satisfies the boundary at a surface of a conical body 
defined by 

'Pc='l'a c +'ZA'l'*c cos 716 (67) 

if the terms p^ are small and quantities of the order of 
pnj 1 or higher are negligible with respect to terms of the order 
of p a . At the surface of the body the velocity must be 
tangent to the surface of the body ; therefore, at each point A 
of the body 

/r*\ dp c _X~ naming 
\»A sm p c dd sin p c 

and, therefore, the boundary conditions can be expressed in 
the approximation considered as 

(5,.)* c =(r, a )* 0c +(^ a )^ S(lW c cosnfl+ 

Sl(®OJ#^(W»coBnfl 

7(S*. f (w0# nsin»J) (YAP. c n sin nd') 
sin p c 

where the subscript p^ indicates quantities at the surface 
of the basic circular body and p c indicates quantities at the 
surface of the conical body considered where the parameters 
(pa), are given by equation (60) and define the shape of the 
shock. Because each term of the right-side summation is 
of the order of (PA 1 , the boundary conditions are 

K)* 0c =o 

’ r (68) 

(MK^).]* 0c =-(-^ 0 ) (iMc 

°c * 

Equations (68) show that the basic flow defined by the 
components r, 0 and is the flow corresponding to a circular 
cone of angle ^'and permits relation of the equation of the 
conical shock to the equation of the conical body. The 
coefficients (tv)., (&„)„, and (w) x for different values of n 
can be determined for a given free-stream Mach number and 
value of (pAc or (PA,’, therefore, from equations (60), 
(65), (67), and (68), the flow field around any conical body 
of the type given by equation (67) can be obtained when the 
terms (PAc are small. 

The determination of the quantities (d,) x , (»*)», and (u?)» 
as functions of p can be obtained from the following equation 



Firubi 11.— Conical coordinate system. 
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(see reference 3) : 


bw 


/. ttA wv H /bw. bv n \ 

sini£d0\ a 2 / a 2 \bip ~ f 'sin ip bd) • 

Equation (69), because of expressions (65), gives 
<=’-(»«); 


(89) 


a 6 f J 


, , , / 2 ”»o , 7-1 ,V'o+ s *o c ° 

cot *+U° +-s- 

ag 

v /,+rr 1 %+N cofc i r v 

«o 4 sin \p 


(70) 


Because of the approximation considered, the entropy 
remains constant in each meridian plane outside of a vorticose 
layer of infinitesimal thickness around the body (see reference 
3). Indeed, 


bS 


bw 


bv r 


by. 


• . VSK/ U U T < 

yR W -1 *’* 9m * 


and 


•v r w sin ip-\-v n w cos i p 
(71) 


bS 

bf 


v n sin ip — =— w 


bS 

be 

bS 


(72) 


From equation (71) there results of the order of 

(W,; therefore, where v*^0, ™ is of the order of ipg and 
can be neglected. Then 
a 2 bS 


7 it bif 


■ 0 


bi) r dw . w bv n , , , . 

= -^bif/- w bvt+sii^ *e +v ' v '- w cot * (73) 


or, in the approximation considered, 

by r 


»»=■ 


bif/ 


(74) 


bv* . 


Then y r + is the radius of the hodograph diagram in 

the plane 6 = Constant (see reference 4) and ( v n )» and 
(v,) n can be obtained from a step-by-step calculation from 
4 r =4 , o/ l to 'f'—foc by means of the equations 


I(y»)it] +- a* = [( t’n),] * cos ( — Aip) + [I?* — (y,)n]+ sin ( — Ap) 


?,)j 


[(p r ) 1 ,]^-.Af=[(y,)#]*sin(— A^)-[^„-(y r ) B ] # cos(— Ap)+(R n ) i ) 

(75) 


where R„ at station ip is obtained from equation (70) and 

[(«).!.-[(»,). +%^] ( 


(70) 


can be calculated from the values of (y r )„, (c«)», and (w), 
at ip. 

The value of [(w)„]*_a* can be obtained from 

K«W» (77) 

-®x can be obtained from equation (71) where 

5'=*S'o+5'i2j^* cos n & (78) 

and Si is independent of ip and can bo determined from 
the equations of the shock from the expression 


S, ~(HX 


°5 


where yj/ s is the inclination of the shock. Then 

Si=p. 0 simA^f+y ro (yr)»+y, 0 (e l .)»+ 

Vr 0 w n sin f +y« 0 w» cos ip (79) 

while, if the quantity Si is neglected and the flow is con- 
sidered potential flow, 


w, 


M* 

1 sin ip 


( 80 ) 


The method presented has been applied to the determina- 
tion, for the condition of zero angle of attack, of the flow 
field around a cone having an elliptical cross section with 
axes in the ratio of about 1 to 3, for which experimental 
data were available at M=1.8, and for a cone of elliptical 
cross section of ratio 1 to 1.88. The calculations have 
been performed in the following way: The bodies are 
shown in figures 12 (a) and 12 (b). The value of ip at 
9=0° is equal to 6.3°, while the value of ip at 6= 90° is 
equal to 18.4-. 

The angle p can be expressed as 


Pc=Po c J rP\ c cos 2 9+ip lc cos 40+^, cos 60+ 


iptc cos 80 + tic cos 100 


(81) 


but 


tan ip c = tan p^- 


S tan p, c 

cosVq,, 


therefore, by determining the value of ip at six points, the 


i 
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following values have been obtained: If ^=10°, 

^=-5.07° 

*t=I.94° 

^=-0.84° 


^=0.43° 


^=-0.16° 

Therefore, a 10° cone at zero angle of attack is assumed as 
the basic body. Reference 5 gives tabulated values for a 
10° cone at 4fi= 1.816 and the calculations have been per- 
formed at this Mach number. The table gives 

^=34.45° 

The entropy variations St are small and are neglected. 
Then, 

<V r )„ s = -0.0145 

(»^=0.0255 

tr,)„ s =1.28o 



r 

i" 

L 



(a) 



< »i ZLliptlcal cane with, exle ratio of 1 to I. 
(b) Elliptical cone wtti«jds ratio of 1 to 128. 
Fiotrai 12.— The elliptical cooes analyzed. 




The values of (s r )„ (p») n , and (io) x between ^ and ^ have 
been determined by means of equations (70), (75), and (80). 
Then ^ , 3 , & s , and ^s 3 have been determined from equations 
(68) where 


(^fL— 2 (".)v=- 2XOj 


6000 


The values of fa obtained are 

fa a =-0.24° 


^=5.1°X10- 3 
fa a = — 1.7°X10~* 

^3=7.4° X lO-o 
lt % =-3 o X10- T 

As is shown from the analysis of the values of and fa, the 
shock is very close to a circular shock wave even for large 
departure of the body from the circular cross section, and 
the effect of the terms corresponding to »=6, 8, and 10 is 
very small. The velocity components at the surface of the 
body are obtained from equations (65) at rp=fa and the pres- 
sure distribution presented in figure 13 is obtained. In the 
same figure, the pressure distributions obtained by using the 
same calculated values of (r r )„ (o*)„ and (w) x for Af=1.81 
around an elliptical cone with a cross section having axes in 
the ratio of 1 to 1.88 are also shown. The conical body 
having an ellipse of axis ratio 1 to 3 has the same cross- 
sectional area as a circular cone of fa— 11°. Its pressure 
drag obtained from this calculation is Cd =0.099 in compari- 
son with 0.12 for the circular cone. The conical body having 
as cross section an ellipse with axes in the ratio of 1 to 1.8S 
has a drag coefficient of 0.103, while the equivalent circular 
cone of ^o=10° 30 f has <7.0=0.115. Therefore, those calcu- 
lations indicate that conical bodies of circular cross section 
have larger drag than cones of elliptical cross section. 



Figum IS. — Pressure distribution around the conical bodies at Af —UU 
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The results obtained agree well with the experimental re- 
sults, also, if the body shape chosen requires large values for 
the angle With the same flow fields (tv)*, (■£>„)„ and 

(«;)„, any other conical shape having two planes of symmetry 
represented by equation (81) when ipo is 10° can be obtained 
at ilf= 1.816. It the flow has only one plane of symmetry, 
only the terms in cos 9, cos 3 9, and so forth, must be con- 
sidered; whereas, if no sy mm etry exists, terms in sin rid and 
cos 7 id must be considered. 

The flow fields defined by (t> r )„, and (w)„ can be ob- 
tained and given in tabulated form without a large amount of 
numerical work as for circular cones, and, therefore, the de- 
termination of conical bodies can be performed without 
difficulty in a very short time. 

When the shock shape is somewhat different from a cone 
having circular cross section, the basic flow field must be 
different from the axially symmetric. However, if the basic 
conical-flow components are expressed in the form 


Equation (17) expressed along the characteristic line of the 
first family X l =tan (jSo+ 930 ) in the plane 2 =Constant is 


1 dV 1 
F 0 dx 


ton a ^ ^ _t_ sin* j3 0 dSi , R , Fi« 
tan /3 °_ + -_ i? 


dw\ 1 tan ftp sin ftp 
dz Vq cos(^o+^ Q ) 


(82) 


while along the second characteristic line Xz=tan(^— fl 0 ) llut 
following equation is valid: 


1 dVt 
V 0 dx 


tonfl «,F| r 


dw 1 1 tan go sin ftp 
— dz V 0 cos (%6to— j8 0 ) 


where B h B 2 , C u and C 2 are defined by equations (20). 
Along each streamline s 0 , 


Vr=v Ta +v rb f(0) 

v n =v Xa +v n J(B) 

w=wj' , (9) 

the basic flow can still be obtained by solving numerically 
the equations of motion in two meridian planes, and, there- 
fore, the basic flow can be determined exactly. For conical 
flow the linearized method can, then, have wide application 
to any form of boundary conditions. 


. dwi &Fi , dSi a 0 2 , S4 v 

dz f dz Wi'o f 

Equations (82), (83), and (84) permit the determination of 
the flow' field by relatively simple procedures. 

Consider, for example, a wing having twist, variable; profile 
distribution, and variable chord, as shown in figure 14. The 
wing can be analyzed by means of the linearized character- 
istics method in the following way: First, the root and tip 
profiles are considered. Section a and section b have differ- 
ent relative thicknesses and chords. 


FLOW FIELDS ABOUND SLENDER BODIES WITHOUT SYMMETRY 

When the conical flow is determined, the method of charac- 
teristics can be applied to the determination of slender 
bodies. The equations used are similar to the equations for 
circular bodies at angles of attack and can be directly de- 
rived from those equations (reference 2). 

Few' values of n are required for the determination of the 
flow' field, and one set of calculations can be used for several 
bodies having the same basic body; therefore, the method can 
be of interest for practical applications. 

QUASI-TWO-DIMENSIONAL FLOW FIELDS 

In many general three-dimensional flow fields of practical 
interest the flow is not too different from a two-dimensional 
flow', and, therefore, the velocity field and entropy field can be 
expressed as in equations (1) and (2) with good practical 
approximation. Flow fields of this kind are found, for ex- 
ample, in wings having plan forms which can be considered 
dose to the two-dimensional type wdtli some tw'ist or a varia- 
tion of thickness distribution along the span. Flow fields of 
this kind can be considered also in some problems in which 
interference between a w'ing and a two-dimensional tail 
(downwash effects) or between a two-dimensional wing and a 
body is considered. In all these problems of practical in- 
terest for the airplane design, the component w in the direc- 
tion of the span of the w'ing can be considered small; therefore, 
equations (1) can bo used and the components iq and v x de- 
pending on the three-dimensional effect can also be considered 
small. 



(a) Wing plan form. 

(b) Wing cross section at a piano j- Constant. 

(c) 8boc]c-wave calculation at the leading edge. 

Fimiua 14.— Determination of the flow aronnd a three-dimensional su person lo wing 
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If the variation from a to b is linear, the properties of a 
two-dimensional cross section at any station e can be obtained 
by means of linear interpolation between the corresponding 
values at a and b. 

The profiles a and b are analyzed by means of two- 
dimensional-flow theory and the characteristic net, and the 
values of the coefficients B u B it C\, and C± are determined 
from two-dimensional considerations. If entropy effects are 
neglected or incorporated in the linearized flow, the coeffi- 
cients B u B. 2i Ci, and C 2 can be determined as for the case 
of two-dimensional potential flow at each point of the axis 
( equations (29). (30) , and (34)) and are constant along char- 
acteristic lines of the first family. Then the linearized flow 
ls defined as the flow that considers the three-dimensional 
effects and the entropy distribution. Therefore, 


u=u n {x,y,z) -\-aUi{x,y,z) j 
!'=Vo{x t y,z) -\-avi(x,y,s) > 
w=awi(x,y,z) ) 


( 85 ) 


where and r 0 are the potential-flow solutions in the plane 
^—Constant and satisfy the boundary conditions in the 
plane z— Constant, u, and r, are the components due to the 
presence of «?, (and of the variation of entropy) , and a can 
he a coefficient, for example, proportional to the. twist dis- 
tribution or to the thickness variation. Because ir 0 and r 0 
are functions .of x, y, and z, while for the basic flow they have 

da Q 
dz 


been determined from two-dimensional considerations, 


( 86 ) 


and are not zero; therefore, equation (22) becomes 

dw__du? u.dw v u , v , dF 
da„~ dz V + d y T = ~T v+ V^ + dz 
or 

dtci ^dl'p . ar t /d Sg . aSA aj 

da, dz dz V dz ^ dz ) yRV a 

dFo . 

where is the variation of the velocity component for the 
basic flow. 

When w L is considered small, in all the flow field the terms 

and ~ can still be neglected in the differential 
Hi oz at oz 

equations along the characteristic lines, and, therefore, equa- 
tions (82) and (83) are still valid. 

At each plane z=Constant, the characteristic net is known; 
therefore, the intersections of the shock wave for the total 
flow with a plane j/=Constant can be determined from char- 
acteristic calculations. If O is a point at the leading edge 
of the wing (fig. 14 (b) j, the'shock wave at O can be obtained 
from shock-wave considerations and from the boundary con- 
ditions because 5 at 0 is known, and at O the shock is two- 
dimensional. Therefore, the velocity components a, and iq 

at O are zero, while — is given by equation (86). If the ’ 

0£n 

plane x—x^ is assumed to be close to the plane r=x 0 , the 
characteristic lines BA and CA can be drawn in any meridian 
plane considered for the basic flow. 


At the point 0, F, is zero and, in the neighborhood of 0 
along the shock L, the velocity can be expressed as 


cro.-ffiii 


(87) 


Now, along the shock wave the direction of the velocity be- 
hind the shock is related to the intensity from the equations 

of the shock wave; therefore, the direction along 

the shock and the shape of the shock as a function of 

are determined when along the shock is known. 

If the velocity F, at a point A of the body is 




( 88 ) 


then the velocity at C and B can be determined as a funct'on 

of from equations (82) and (83), because the value of uq 

at- A is given from equation 86 and is known, and the value 

of can be obtained from the value of w, at A in several 
dz 

planes z=Constant. 

If in equations (82) and (83) the values of F,, Si, and p, 
at B and C are expressed by means of equation (87) and of the 
equations of the shock waves which give the coefficients of 
the expressions 


dF, 

a VidL 

dx 

= a l dx 

d‘Pi i 

a<p, ar, 

dx 

ar, a r 

as, 

as, ar, 


dx dF, dx 


then equations (82) and (S3) give two relations between 

and and, therefore, and 4V can be determined. 
ds dL ’ ’ ds dL 

The equation of the shock wave can relate u as a function 

of v, or Vi as a function of <pi without the necessity of the 

component tc, because the component w is proportional to the 

inclination y of the tangent to the shock with the plane 

x= Constant, and Ft and p t ar e functions of Mi cos ij, but 

Mi cos y^=A/j in the approximation considered 

here (fig. 14 (c)). The components u and v at B and the 
position of B having been determined, in each meridian 
plane z= Constant, the intersection of the shock in 
the plane ;/= Constant is obtained and tr at B is determined. 
Then a point D is interpolated in each meridian plane and 
the point E is obtained. Then point F is determined. In 
order to obtain (■ w ) r , the streamline DF' for the basic flow 

must be drawn and interpolated between F' and E. 

In a similar way, all the flow field can be obtained. The line 
TT' defining the plan form in figure 14 (a) must be outside of 
the Mach conoid from T. Because few points along each 
profile are required, the largest amount of work for such a 
calculation is represented by the construction of a basic 
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characteristic net which permits obtaining points which 
simplify the determination of’tr. By changing the value of 
the coefficient a, different thickness distributions or different 
twists can be considered. The new distributions must be 
obtainod by changing proportionately the variation of thick- 
ness or twist with respect to the basic w r ing and by varying in 
proportion the value of a. 

CONCLUDING REMARKS 

The method of characteristics for supersonic flow has been 
simplified by assuming that one of the velocity components 
or the effect on the velocity components due to variation 
•of one physical parmater is small, so that the square of the 
velocity components considered small can be neglected. 
By means of this simplification, the flow field can be repre- 
sented as the superposition on a basic flow field (which is not 
linear and must be determined by the method of character- 
istics) of linearized flow fields which are defined by a differ- 
ential equation with variable, but known, coefficients. 

The calculations of these linearized flow fields can be per- 
formed along the characteristic net of the basic flow field. 
The method has been applied (a) to the two-dimensional flow 
with entropy gradient, which has been transformed to a basic 
potential flow on which a linearized flow due to the entropy 
gradient is superposed, (b) to axially symmetric problems 
where conical or cylindrical flows are considered as the 
basic flow, (c) to the determination of the flow field around 
cones or slender bodies without axial symmetry, and (d) to 
particular three-dimensional flows which can be simulated 
as a basic two-dimensional flow on which three-dimensional 
linearized flows are superposed. Application (b) permits 
obtaining in a simple way the flow field around bodies of 
revolution without using linearized theory and indicates the [ 


possibility of using tabulated values for such determinations. 
Application (c) permits the determination of flow' fields 
not yet determined by the method of characteristics. Any 
such conical flow's can be determined by using tabulated 
values that can be obtained as for cones of circular cross 
section at small angles of attack. The application in (d) 
can bo of interest for wings of approximately tw'o-dimensional 
form having tw'ist or thickness variation along the span and 
for interference problems. 

Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Ya., July 24, 1951. 

REFERENCES 

1 . Sauer, Robert: Supersonic Flow About Projectile Heads of Arbitrary 

Shape at Small Incidence. R..T.P. Translation No. 1573, 
British Ministry of Aircraft Production. (From Luftfalirt- 
forschung, vol. 19, no. 4, May 1942, pp. 148-162.) 

2. Ferri, Antonio: The Method of Characteristics for the Determina- 

tion of Supersonic Flow Over Bodies of Revolution at Small 
Angles of Attack. NACA Rep. 1044, 1951. (Supersedes 
NACA TN 1809). 

3. Ferri, Antonio: Supersonic Flow Around Circular Cones at Angles of 

Attack. NACA Rep. 1045, 1951. (Supersedes NACA TN 
2236.) 

4. Ferri, Antonio: Elements of Aerodynamics of Supersonic Flows. 

The Macmillan Co., 1949. 

5. Staff of the Computing Seotion, Center of Analysis (Under Direction 

of Zdenek Kopal): Tables of Supersonic Flow Around Cones. 
Tech. Rep. No. 1, M. L T., 1947. 

6. Busemann, A.: Conical Supersonic Flow With Axial Symmetry. 

R.T.P. Translation No. 1698, British Ministry of Aircraft 
Production. (From Luftfahrtforschung, vol. 19, no. 4, June 5, 
1942, pp. 137-144.) 



